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Abstract 


An exact time-dependent solution for the wave function ip(r, t) of a 
particle moving in the presence of an asymmetric rectangular well/barrier 
potential varying in one dimension is obtained by applying a novel for this 
problem approach using multiple scattering theory (MST) for the calcula¬ 
tion of the space-time propagator. This approach, based on the localized 
at the potential jumps effective potentials responsible for transmission 
through and reflection from the considered rectangular potential, enables 
considering these processes from a particle (rather than a wave) point of 
view. The solution describes these quantum phenomena as a function of 
time and is related to the fundamental issues (such as measuring time) of 
quantum mechanics. It is presented in terms of integrals of elementary 
functions and is a sum of the forward- and backward-moving components 
of the wave packet. The relative contribution of these components and 
their interference as well as of the potential asymmetry to the probability 
density \ip(x,t)\ 2 and particle dwell time is considered and numerically 
visualized for narrow and broad energy (momentum) distributions of the 
initial Gaussian wave packet. The obtained solution is also related to 
the kinetic theory of nanostructures due to the fact that the considered 
potential can model the spin-dependent potential profile of the magnetic 
multilayers used in spintronics devices. 
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1 Introduction 

The time-dependent aspects of reflection from and transmission through a po¬ 
tential step/barrier/well raised several questions that have not yet been com¬ 
pletely clarified and have recently acquired relevance in view of renewed interest 
in the fundamental problem of measuring time in quantum mechanics (see HD- 
The tunneling time can serve as an example (see, e.g. a review 0 ). 

The mentioned phenomena are less surprising when we think of a wave be¬ 
ing, e.g., reflected from a downward potential step. In the stationary case, these 
quantum phenomena easily follow from standard textbook analysis, which re¬ 
duces to solving the stationary Schrodinger equation by matching the wave 
function of a plane wave of energy E and its derivative across the potential 
jumps. However, in this case, there are no real transport phenomena, i.e. in the 
absence of energy dispersion, A E = 0, the particle time of transmission through 
or of arrival (TOA) to the potential jumps is indefinite (At ^ H/AE). 

These processes are more surprising from the particle point of view, and it 
is interesting to verify the mentioned non-classical phenomena by considering 
the time-dependent picture in a realistic situation, when a particle, originally 
localized outside the potential wcll/barrier, moves towards the potential and 
experiences scattering at the potential jumps. In order to describe these time- 
dependent processes, the corresponding time-dependent Schrodinger equation 
with a rectangular potential should be solved, which is much more involved 
compared to the conventional stationary case. From the particle point of view, 
it also seems desirable to be able to apply the multiple scattering theory (MST) 
to this principally important situation. The MST is conventionally formulated 
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in terms of a ’’free” particle Green function (propagator) modified by the scat¬ 
terings at the inhomogeneities such as the interfaces between different media. It 
is clear that the scattering events occur in the interface area and are stipulated 
by some fairly localized potentials. Thus one faces an interesting problem of 
finding a localized potential responsible for particle scattering from a potential 
inhomogeneity. 

In addition, there is one striking and classically forbidden counterintuitive 
(and often overlooked) effect even in the simplest ID time-dependent scattering 
by the mentioned potentials. A wave packet representing an ensemble of parti¬ 
cles, confined initially (at t = to), say, somewhere to the region x < 0, consists 
of both positive and negative momentum components due to the fact that a 
particle cannot be completely localized at x < 0 if the wave packet contains 
only p > 0 components. One would expect that only particles with positive 
momenta p may arrive at positive positions x > 0 at t > to. However, the wave 
packet’s negative momentum components (restricted to a half line in momen¬ 
tum space) are necessarily different from zero in the whole x space (—oo -r- oo), 
represent the particles’ presence at x > 0 at initial moment of time to, and, 
therefore, may contribute, for example, to the distribution of the particles’ time 
of arrival (TOA) to x > 0 [3J H]. It is worth noting that the contribution 
of the backward-moving (negative momentum) components in the initial-value 
problem is in some sense equivalent to the contribution of the negative energy 
(evanescent) components in the source solution [3]. Thus, the correct treatment 
of some aspects of the kinetics of the wave packet (even in the ID case and even 
for ’’free” motion) becomes a nontrivial problem and is closely related to the 
fundamental problem of measuring time in quantum mechanics, such as TOA, 
the dwell time, and tunneling time. 

On the other hand, the mentioned rather academic (but fundamentally im¬ 
portant) problems have acquired reality and significance due to important prac¬ 
tical applications in the newly emerged field of nanoscience and nanotechnol¬ 
ogy. Rectangular potential barriers/wells may often satisfactorily approximate 
the one-dimensional potential profiles in layered magnetic nanostructures (with 
sharp interfaces). In such nanostructures, the giant magnetoresistance (GMR) 
[5] and tunneling magnetoresistance (TMR) [5] effects occur. These effects, 
which stem, particularly, from quantum mechanical spin-dependent electrons 
tunneling through potential barriers or their reflection from potential wells, 
have led to very important commercial applications of spintronic devices. 

The solution to the time-dependent Sclrrodinger equation can be obtained 
with the help of the spacetime propagator (Green’s function), which has been 
conveniently calculated by the path-integral method. The list of exact solutions 
for this propagator is very short. For example, there is an exact solution to 
the space-time propagator by the path-integral method in the one-dinrensional 
square barrier case obtained in [7], but this solution is very complicated, implicit 
and not easy to analyze (see also [51 1§1 !T0]). 

Recently, we have suggested a method m for the calculation of the space- 
time propagator which is based on the energy integration of the spectral density 
matrix (discontinuity of the energy-dependent Green function across the real en- 
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ergy axis). The energy-dependent Green function is then easily obtained for the 
step/barrier/well potentials with multiple-scattering theory (MST) using the 
effective energy-dependent potentials found in m, which are responsible for re¬ 
flection from and transmission through the potential step. The obtained 5-like 
potentials describing the quantum scattering from a potential step provide a 
clear picture of the particle’s scattering taking place at the interfaces and make 
it convenient to calculate the energy-dependent Green function and the space- 
time propagator especially when the scattering from more than one interface 
needs to be accounted for and there are other sources of scattering by the point¬ 
like scatterers. Such a situation is typical for real multilayers with disordered 
interfaces mmm and for the Casimir effect m- An important advantage 
of our approach to propagator calculation is also that it allows for a natural 
decomposition of the general initial wave function evolution in time into both 
the forward- (p = hk > 0) and backward-moving (negative momentum p < 0) 
components and an analysis of the contribution of both these terms and their 
interference to particle reflection and transmission. This approach has been fur¬ 
ther applied to the analysis of the time-dependent properties of the scattering 
by the imaginary step JT6] (related to calculation of the particle time of arrival) 
and rectangular symmetric barrier/wcll potentials pT] ITS] . 

In this paper, we generalize our approach to the solution of the time-dependent 
Schrodinger equation in the case when a particle moves towards a rectangular 
asymmetric (spin-dependent) potential. Tunneling through asymmetric poten¬ 
tials has already important applications in semiconductor heterostructures. The 
asymmetric (spin-dependent) rectangular potential barrier/well can also model 
the potential profile of the magnetic tlrreelayer switched from the parallel con¬ 
figuration of magnetic layers (symmetric potential) to the anti-parallel config¬ 
uration of layers. Although this potential, which models the spin-dependent 
potential profile in magnetic nanostructures, changes only in the x (perpendicu¬ 
lar to interfaces) direction, the system under consideration in this paper is a real 
three-dimensional one. A simple exact solution for the time-dependent propa¬ 
gator in terms of integrals of elementary functions is obtained, which is valid for 
both the well and barrier cases. This solution fully resolves the corresponding 
time-dependent Schrodinger equation and provides exact analytical expressions 
for the wave function ip{r,t) in the spatial regions before, inside and after the 
potential with account for the backward-moving terms caused by the negative 
momentum components of the initial wave function. It is important that the 
obtained solution allows for numerical visualization of the observables defined 
by the wave function if>(r,t) in the mentioned spatial regions. Thus, the corre¬ 
sponding probability densities |^(r, t)\ 2 are analyzed and numerically visualized 
for the Gaussian initial wave packet with special attention to the counterintu¬ 
itive contribution (see M) of the backward-moving wave packet components 
and the potential asymmetry. We show (and visualize) that the contribution 
of the backward-moving components of the wave packet is small in the quasi- 
classical case but is otherwise important. It is also shown that the influence 
of the potential asymmetry is more pronounced when the contribution of the 
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backward-moving wave packet components is essential. The dwell time, which 
characterizes the average time spent by a particle in the potential region and 
is related to the enduring quantum physics problem of calculating the tunnel¬ 
ing time, is considered for the asymmetric rectangular potential. The obtained 
results can also provide a foundation for a kinetic theory of nanostructures. 

2 Multiple-scattering calculation of the space- 
time propagator and time-dependent solution 
for the Schrodinger equation 

We consider a particle moving toward the following asymmetric one-dimensional 
rectangular potential of the width d placed in the interval (0 < x < d) 

V(x) = [O(x)-O(x-d)]U + 0(x-d) A, (1) 

where 9{x) is the Heaviside step function, and the potential parameter U can 
acquire positive (barrier) as well as negative (well) values. The wave packet, 
modeling a particle, will approach a potential m from the left (where the parti¬ 
cle potential energy is zero) and the parameter A is supposed to be non-negative 
(A > 0). With the potential dU) we can model, e.g, the spin-dependent poten¬ 
tial of a threelayer, which consists of a spacer (metallic or insulator) sandwiched 
between two magnetic (infinite) layers. An asymmetry (spin-dependence) of the 
potential (CD is defined by the parameter A via the electron spectrum in different 
magnetic layers as 

*<(-£; k ||) = fc(-E;k||),fc>(-E:k||) = fc A (E;k||), 

fc (^ k ll) = \Jw E ~ ^’^kN) = ]/|?(^- A )- k ||> 

(E; k M ) = k<(E; k,,) = k u (E- k„), k u (E; k,,) = ^ (E - U) - kjj, (2) 

where ( < )(-E'; k||) and k||) are the perpendicular-to-interfaces (lo¬ 

cated at x = 0 and x = d) components of the particle wave vector k to the right 
(>) or to the left (<) of the corresponding interface, while k|| is the parallel- 
to-interfaces component of an electron wave vector, which is conserved for the 
sharp interfaces under consideration. The two-dimensional vector k|| defines 
the angle of electron incidence at the interface. 

From the particle propagation point of view, the partial reflection from and 
transmission through a potential inhomogeneity may be explained by the quan¬ 
tum mechanical rules of computing the probabilities of different events. These 
rules represent the quantum mechanical generalization of the Huygens-Fresnel 
principle and were introduced by Feynman as the path-integral formalism [IS]. 
It states that a wave function of a single particle moving in a perturbing poten- 
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tial V (r, t) may be presented as 


^(r,i) = j dr'K{r,t\r',t 0 )^(r',t 0 ). (3) 

Equation m shows (in accordance with the Huygens-Fresnel principle) that the 
wave function ^(r, t) at the spacetime point (r, t) is the sum of the contributions 
of all points of space where the wave function ^(r', to) at t = to is nonzero. The 
propagator K( r, t; r', t 0 ) is the probability amplitude for the particle’s transition 
from the initial spacetime point (r',fo) to the final point (r,t) by means of all 
possible paths. It provides the complete information on the particle’s dynamics 
and resolves the corresponding time-dependent Schrodinger equation. 

Thus, the problem is to find the propagator K(r, t; r', t 0 ) for the given po¬ 
tential V (r, t). In some cases, for example when the potential is quadratic in the 
space variable, the kernel K (r, t; r', to) may be calculated exactly. In the case 
when the potential changes smoothly enough, a quasi-classical approximation 
can be employed. It is not, however, the case for the singular potential ©■ 
According to m , the time-dependent retarded (operator) propagator K(t\t') 
6(t — t') exp [~j-H(t — t')] can be calculated with the use of the following defi¬ 
nition 


OO 

K(t-t’) = 0(t-t')^- [ e-i E ( t - t ''>[G(E + i£)-G(E-i£)}dE,e^+ 0, (4) 
27T J 

— OO 


where 


G{E) 


1 

E — H 


(5) 


is the resolvent operator, E stands for the energy and H is the Hamiltonian of the 
system under consideration. Correspondingly, G{E ± *e) = G± (A) defines the 
retarded (G + ) or the advanced ( G ~) Green function. The E-resolving Fourier 
transformation © is useful for the calculation of the propagator K(t;t') when 
the Green functions G ± (E) may be found for each value of E, i.e. when the 
considered processes are energy-conserved as is the case considered in this paper. 

We are looking for the spacetime propagator A'(r,f;r',0) =< r|A'(t;0)|r' > 

, defining the probability amplitude for a particle’s transition from the initial 
point (r',0) to the final destination (r ,t) in the presence of the potential ©. 
For the considered geometry, it is convenient to present the r-representation of 
the Green function with the Hamiltonian H , G(r,r';A) = < r| E ^ H |r' >, as 
follows 

G(r,r/;E) = I^n^G^, a/; £; k M ), (6) 

A kn 


where p = ( y , z) is a two-dimensional parallel-to-interface vector and A is the 
area of the interface. Thus, the problem is reduced to finding the one-dimensional 
Green’s function G(x , a/; E; k||) dependent on the conserved particle energy and 
parallel-to-interface component of the wave vector. In the following calculation 
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of this function we will suppress for simplicity the dependence on the argument 
k| |, which will be recovered at the end of calculation. 

We showed in m that the Hamiltonian corresponding to the energy-conserving 
processes of scattering at potential steps can be presented as 


H = H 0 + Hi(x-,E), 

H i (x-,E) = Y, H i( E )t(x- x s). (7) 

S 


Here, Hi(x;E ) describes the perturbation of the ’’free” particle motion (defined 
by H 0 = localized at the potential steps with coordinates x s (in the 

case of the potential ©, there are two potential steps at x s = 0 and x s = d ) 


H? > (E) = ^[v° > (E)-v s < (E)}, 

H! < (E) = ^[v%(E)-v s > (E)}, 
2 ihvUE^E) 


( 8 ) 


where is the reflection (from the potential step at x = x s , s € {0, d}) 

potential amplitude, the index > (<) indicates the side on which the particle 
approaches the interface at x = x s : right (>) or left (<); H? >< (E) is the 
transmission potential amplitude, and the velocities = hk s > ^{E)/m, 

(k>(<){E) are given by ©). Note that the perturbation Hamiltonian Hf ’s 
dependence on k|| (which is omitted for brevity) comes from Eq. ©. 

The perturbation expansion for the retarded Green function G + (x,x'’,E) in 
the case of the rectangular potential m, which can be effectively represented 
by the two-step effective Hamiltonian 0, reads for different source (given by 
x ') and destination (determined by x) areas of interest as (see also (17, IS) 


G+(x,x'-,E) 
G+(x,x'-E) 
G+{x,x'-,E) 
G+(x,x'-E) 
G+(x,x'-E) 


G^(x,d;E)T + {E)G^{0,x';E),x' < 0 ,x>d, 

G+(x,0;E)T + (E)G+{d,x';E),x' > d,x < 0, 

G$(x, 0; E)T' + {E)Gq (0, x'; E) + G$(x, d; E)R'+(E)G£{ 0, x';E),x / < 0,0 < x < d, 
G^{x,0-,E)T' + (E)G^(0 1 x';E)+ G^(x,0;E)R'+(E)G^{d,x';E),0 < x 1 < d,x < 0, 
G+{x,x'-E) + G+(x,0;E)R+(E)G+(0,x';E),x' < 0,x < 0, 

(9) 
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where the transmission and reflection matrices are 


T+(p , _ r^+(E)G+(rf,0;E)r°+(E;) 

( ' D+(E) 

T ' + (E) = ^^y’ R ' + ( E ) = T* + (E)G+(d,0;E)T'+(E), 

R+(E) = T°+(E) + T >< {E)G ° ( °’ d; E ^ T < + W G t K °: E ) T >< W 

D+(E) = 1 - G+(d, 0; £)T°+(£)G+(0, d; E)T*+(E). (10) 

The one-dimensional retarded Green function Gq(x, x'\ E) corresponding to a 
free particle moving in constant potential V(x) =0 or V(x) = U{ or A) is (see, 

e-g. ED]) 

Gq (x, x'\ E) = ih 2 ^ E j exp[ik(E)\x - x'\],V{x) =0, 

Gq(x,x';E) = ™ ex Pl ik u(A){E)\x - x'\],V{x) = U{ or A), (11) 

where the wave numbers are determined by 0. The scattering (at the step 
located at x = x s ) t-matrices are defined by the following perturbation expan¬ 
sion: 


T S {E) = H?(E) + H?(E)G 0 (x a ,x a -,E)H!(E) + ... 
H?(E) 

1 — G 0 (i s , x s ; E)H? (E) ’ 


( 12 ) 


where H-(E) and the interface Green function Go(x s , x s ; E) are defined differ¬ 
ently for reflection and transmission processes El: the step-localized effective 
potential is given by Eq. © and the retarded Green functions at the interface 
for the considered reflection and transmission processes are, correspondingly, 


G J>«) ( x s, x s',E) = l/ihv a >{<) (E) 

Gj><(i s ,a: s ;£) = 1 /ih^v%(E)v^(E) (13) 

in accordance with m- 

From ©, (1121) and (1 1j f) . we have for the reflection T^ < ^(E) and transmission 

T^(E) t-matrices, used in (flOl) (s G {0, d}), corresponding to the retarded 
Green function and scattering at the interface located at x = x s G {0, d}, 

E >l<)^ E ^ = *^>(<) r >(<)> 

TS+(E)=ihy/^t s , (14) 


where r^^JE) and t s (E) are the standard amplitudes for reflection to the right 










(left) of the potential step at x = x s and transmission through this step 


j US _ L,s b s — b s 

rm = tr J ^>r° < (E)= < > 


*£+*£’ <v 


k s > + k s K ’ 


t s (E) = 


k> + k s K 


( 15 ) 


and the argument E in the wave vectors is omitted for brevity. 
Using Eqs. (El,©, (fTUl) . (fill) , (fill) and ESI, we obtain 

G+{x,x'-E) = o m r _ e ik ^ x - d h(E)e- ikx ', x' < 0 ,x>d, 

ih 2 ^/kkA 


G + (x,x';E) = 


ih 2 \fkk\ 


e~ ikx t{E)e ik ^ x '- d \x' > d,x < 0, 


G + (x, x';E) = - 


G+{x,x'-,E) = 
G+(x,x'-,E) = 


iH 2 y/kk u - 
m 

ih 2 yjkk u - 


e ikuX t'(E)e~ ikx ' + e ~ ikuX r' (E)e~ ikx> 


G'(E)t 


—ikx^J ( tp\ „—ik. 


(E)e 


m 


ik\x—x'\ 


+ r{E)e- ik ( x+x) ] ,x < 0,x’ < 0, 


ih 2 k L 

where the transmission and reflection amplitudes are defined as 

^kk A k u e lkud , 2y/kk^(k A + k u ) 

t(E) = - —— - ,t (E) = 


, x 1 < 0,0 < x < d, 
, a; < 0, 0 < a :' < d, 

(16) 


d(E) 


d(E) 


, 2^kk^(k u - k A )e 2lkud (fc - k u )(k A + k u ) - (fc + k u )(k A - k u )e 2zkud 

r \ E ) = - J7TT\ - H E ) = - 


d{E) 

d{E) = (fc + k u )(k A + k u ) - (.k - k u ){k A - k u )e 


d(E) 


2 ik u d 


(17) 


We remind that k, k u and k A are the perpendicular-to-interface components of 
the particle wave vector in different spatial areas which depend on the energy 
E and k|| as indicated in m, and k|| is the parallel-to-interface component 
of this vector which is conserved for the considered specular scattering at the 
interfaces. Using the same approach, it is not difficult to obtain the Green 
function G + (:r, x'\ E) for other areas of arguments x and x'. 

The transmission probability \t(E)\ through and reflection probability \r(E)\ 
from the asymmetric potential ED, which follow from E3 for real k u and k A , 
are given by 


\t(E)\ 2 
HE) | 2 


4kk 2 k A 

( k + k A ) 2 k 2 + ( k 2 - k 2 )(k\ - fc 2 ) sin 2 {k u d) ’ 
k 2 {k - k A ) 2 + (fc 2 - k 2 )(k A — fc 2 ) sin 2 (k u d) 
(k + k A ) 2 k 2 + (fc 2 — k 2 )(k\ — fc 2 ) sin 2 (fc„d) 


(18) 


Note that when k u d = mr (n is integer), the resonance transmission (\t(E)\ 2 = 1 
and \r(E)\ 2 = 0) happens only for a symmetric rectangular potential with A = 0 
{k A = k). 
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In accordance with the obtained results for Green’s functions, we will con¬ 
sider the situation when a particle, given originally by a wave packet localized 
to the left of the potential area, i.e. at x' < 0, moves towards the potential ©■ 
We also choose A ^ 0, which corresponds to the case when, e.g., the spin-up 
electrons of the left magnetic layer (a/ < 0) move through the nonmagnetic 
spacer to the right magnetic layer (x > d) aligned either in parallel (A = 0) or 
antiparallel (A > 0) to the left magnetic layer. At the same time, the amplitude 
U in the potential © may acquire both positive (barrier) and negative (well) 
values. 

From Eqs. m we see that G + {x,x'\E) = G + (x',x-,E), and, therefore, the 
advanced Green function G~ (x,x' \E) = [G + (x', x-, E)]* = [G + (x, x'; E)}* (see, 
e.g. 123). Thus, the transmission amplitude © is determined by the imaginary 
part of the Green function and can be written as 


OO 

K(x, t;x',to) = — 6(t — to)— J dEe~'K E ^ t ~ to ^ lmG + (x,x'-,E). 


(19) 


Formulas (USD - © present the exact solution for the particle propagator 
in the presence of the potential © in terms of integrals of elementary functions 
for a given angle (k| |) of a particle’s arrival at the potential ©. Thus the Green 
function and propagator are dependent on the additional argument k||, i.e. ac¬ 
tually we have obtained the solution for G + (x , x'; E\ k||) and K(x, t; x ’, to; k| |)- 
It should be kept in mind that the wave numbers © and, therefore, the quan¬ 
tities t(£';k||), t , (£';k||), r'(E: k||), and r(i?;k||) in (fTTll are different in the 

fi2 kf,/2m oo 

J dE and J dE energy integration areas: in the former case, k(E- : k||) 
— 00 ft 2 k^/2m 

and k/\(E;k\\) (A ^ 0) should be replaced with ik(E; k||) and ikA{E;k^), 
where fc(A;k||) = .^/kj^ — 2mE/h 2 and A’A(F';k||) = ^/k^ + 2m(A — E)/h 2 . 
At the same time, for energies E < fi 2 k^/2m, the wave number k u = ik u , 
k u = ^/kj^ + 2m(U — E)/h 2 , for U >0 (barrier), but for U < 0 it is real, 
i.e. k u = yj2m{E + \U\)/h 2 — k^, if E > H 2 k^/2m — \U\ and k u = ik u , 
k u = ~ 2m(E + \U\)/h 2 if E < h 2 k 2 ^/2m — |t/|. It follows that the 

’’free” Green function Gq(x,x'-,E) = j^e zk \ x ~ x I is real in the energy inter¬ 
val (— oo 4- h 2 k 2 |/2m) and, therefore, does not contribute in this interval to 
the corresponding ’’free” propagator Kq(x, t; x', to) defined by (USD . It is also 
remarkable that for energies E < H 2 ’k. 2 ^/2m the imaginary parts of the Green 
functions vanish in all spatial regions, as is seen from definitions ( 1161 ) and © 
(e.g., Im t(E) = 0 and Imr(Fl) = 0 for E < h 2 k 2 ^/2m). Therefore, the energy 
interval (—oo-^/i 2 k 2 | /2m) does not contribute to the propagation of the particles 
through the potential well/barrier region. 
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From Eqs. ©, © and © we see that the dependence of the Green 
function on E and k|| comes in the combination E — trk^/2TO, and, therefore, 
it is convenient to shift to this new energy variable, which is the perpendicular- 
to-interface component of the total particle energy. Thus, accounting for © 
- © and that for the new energy variable the energy interval (—00 4- 0) does 
not contribute to the propagator, we have for t > to 


K(x,t;x',t 0 ;k||) = 


K(x,t;x’,t 0 ; k||) = 


_ i 1 

g h 2m 




Tlh 


dEe-* E{t ~ to) - 


Vv{E) 


Re 


\/va(E) 


t{E)e ik *( E )( x ~ d ) e ~ ik ( E ) x ' 


x' < 0, x > d, 


a 2 kS 00 

e~h 2 m (*—* 0 ) r i . n 1 

- / dE P -nEjt-to) Re 


7 tH 


\/v(E) I tJv u (E) L 


p —ik(E)x' r 

- t'(E)e iku{E)x + r'(E)e~ iku{E): 


x' < 0,0 < x < d, 


/v(a;,t;a; , ,to;k||) = 


e ft 2m (* *0) /' _ , N 1 


dEe-i E{t - to) —-Rc[e iKE) \ x ~ x '\ + r(E)e- ik< ' E)< ' x+x, ' l },x' < 0,x < 

7r ri J V \E) 

0 

( 20 ) 


where (see ©) 


v{E) = Hk(E)/m,k(E) = i/2 mE/h 2 , 
va(E) = hkA(E)/m , k^iE) = \j2m{E — A )/h 2 , 
v u {E) = hk u (E)/m,k u (E) = \/2m(E- U)/H 2 . (21) 


The transmission and reflection amplitudes t(E), r(E), t'(E ) and r'(E ) in (1201) 
are defined by © with the wave numbers ©. 

It is easy to verify that the integration over E and k|| (according to (|6|l) 
of the first term in the last line of © results in the known formula for the 
space-time propagator for a freely moving particle 


K 0 (r,t-r',t 0 ) = 6(t - t 0 ) 


m 


3/2 


2nih(t — to) 


exp 


im( r — r') 2 
2H(t — to) 


, x < 0, x' < 0. 


( 22 ) 

The obtained results for the particle propagator completely resolve (by means 
of Eq. (j3j ) the time-dependent Schrodinger equation for a particle moving under 
the influence of the rectangular potential ©. 
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3 Time-dependent probability density of finding 
a particle in different spatial regions 

Using Eqs. ©, (USD and OSD. we can present the wave function in different 
spatial regions at t > to as 




(23) 


Here 


OO 


ip<(r,t) = 


y/2nhj 


dEe~^ E{t ~ to) 




t*{E)e- ik ^ E)( - x - d ^ < (E-p,t),x > d, 




.'(^ p iku( E )x U(E)e- ifc “(' B ) 


p,t),0 < x < d, 


ip<(r,t) = 




dEe~i E(t ~ to) 


V<(E) 


'* ( F!\p~ ik u( E ) x -I_ r ■'* ( F,\p ik u( E ) 


[i t'*(E) 


+ r'*{E)e lk ^ E)x ]ip < {E;pJ), 0 < x < d, 


OO 

ip>(r,t) = —;L= [dEeriiW-* o) 1 [e ifc ( £ ) x + r(E)p- ik ( E ) x ]il h (E: p.t ), .r <f 0, 

y/27rhJ Wv(E) 


ip<{r,t) = 


y/2nh 


u 
OO 

f dEe-^ E ^~ to) * _ [e~ zfc(g)x + r ,*( jB ) e ifc(s)x]^<( jB ; p ^ x < o 

J V V (E) 

(24) 


and r = ( x,p ). The wave function in the U-representation ip > t < p ) {E\ p,t) is re¬ 
lated to its /c-representation V ; > ( < ) [fc(-E); p,t)] as 

tp>{E;p,t) = -^===^>[fc(S);p,t],V’<(S;P,t) = ^/2 ttL(U) 

-0>[fc(£;);p,t] = J dp'K(p,t;p’,t 0 ) J dx'e~ lk{E)x ', p',t 0 ), 
ip<[k(E);p,t] = J dp'K(p,t;p' ,to) j dx'e lk< ^ E)x '^{x',p',t 0 ), 

i ■ ?i 2 k 2 

K(p,t;p',t 0 ) = 7 ^ ex Ph^^( f -t 0 )]e lk " {p_p,) = 


_IS_ ex p[— 

2nih(t — to) 2ih(t — to) 

(25) 


where K(p,t\ p',to) is the ’’free” propagator in the parallel-to-interface (y,z) 
plane (see ©, OSD and dZZD ). 
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It can be verified that the wave functions EH) and their derivatives are con¬ 
tinuous at x = 0 and x = d. To be definite, we assume that for positive energies 
k(E) = yj2mE/h 2 > 0 and, therefore, [k(E)\ p, t] is related to the compo¬ 
nent of the initial wave function , p', to) corresponding to propagation to the 
right along the x axis, and, accordingly, if) < [k(E ); p, t] represents propagation to 
the left. When the potential V(x) ^ 0, integration over x' in (E5l) is restricted 
to the negative semispace (x 1 < 0), as it follows from the expressions (l20l) for 
the particle propagator. 

The result, given by Eqs. (l23l) . (l24l) and (j25)l . indicates that, generally, the 
contribution of the wave function, originated at t = to to the left of the potential 
© (a/ < 0), to the wave function in the region of the potential (0 < x < d) and 
to the right of it {x > d) comes at t > t$ from both: the components moving 
to the right, i/j > , and to the left, ip <. This rather paradoxical result follows 
from the fact that if the initial wave packet has the non-negligible negative 
momentum components (restricted to a half line in the momentum space), the 
corresponding spatial wave function is different from zero in the entire x-region 
(—oo, oo), interacting with the potential even at t < to , and is thus modified by 
this interaction (see also EH, 0). As a result, the backward-moving components 
contribute to the behavior of the wave function at t > to in the spatial regions 
to the right of the original wave packet localization. 

Consequently, the probability density of finding a particle in the spacetime 
point (r, t), t.)\ 2 is determined by the forward- and backward-moving terms, 

as well as their interference: 


\4’{r,t)\ 2 = |-0>(r, i)| 2 + |-0c(r, £)| 2 + 2 Re^>(r, t)^<(r, t). (26) 

Equations (l24l) - (l26l) generally resolve the problem of finding a particle in the 
spatial region of interest at time t for a given initial wave function ip(r',to). 
These equations can be used for numerical modeling of the corresponding prob¬ 
ability density in the different space-time regions (see below) and for determining 
some characteristics of the particle dynamics under the influence of the potential 
©• 

In order to estimate the actual contribution of the backward-moving and 
interference terms to the obtained general formulas, we should consider a phys¬ 
ically relevant situation as to the initial wave packet. Let us consider the case 
when the moving particles are associated with a wave packet which is initially 
sufficiently well localized to the left of the potential ©• Thus we now consider 
the problem for a particular case of the initial state corresponding to the wave 
packet 




( 27 T ( 7 2 ) 3 / 4 6XP 


(? - Vi? 

4<t 2 



,Xi < 0,fc* > 0, 


(27) 


located in the vicinity of r,; = ( Xi, pi ) and moving in the positive x direction 
with the average momentum p, = Ek. ;, ki = k x > 0 (k,, = (fc*, k®,), r 7 = (x',p')). 
Thus, we consider a general situation, when a particle, associated with the 
wave packet m, comes to the potential m from the left with the positive 
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perpendicular-to-interface momentum component Hki > 0 at the angle defined 
by the parallel-to-inteface momentum component /ikj . Now, we can perform 
integration over spatial variables x', p\ as it follows from (l25l) and (l27l) .The 
result is 


*P>(E) 

ip<(E) 

C\\(p,t) 


C\\(p,t)ip > (E),'ijj < (E-, p,t) = C\\(p,t)ij} < (E), 

1 p m 2 ] l/4,i[t i -t(E)li v -|t.-t(S)] 2 » 2 

tHv(E) 

(2na 2 ) 1/ ' l e i[ki+kiE ' >]xi e~ lki+k{E ' )]2a2 , 


\Jnhv(E) 

<2 


7 r ih(t — to) + 2ma 2 


exp[- 


(p- pi - 2«k[|C t 2 ) 2 ?7 

2 ih(t — to) + Ama 2 


-]e ik i l f’ i e”( k l l <T ) 2 


(28) 


where the factor C\ \ ( p,t) defines the dependence on the parallel-to-interface 
components of the vectors involved. Thus, the forward- and backward-moving 
components of the wave function (r,t) CH) for the initial wave packet (1271) 
reduce to the one-dimensional integral over energy E with the energy-dependent 
functions ip > i < ' ) (E) and the common factor C\\(p,t). 

We note that 

J dp\C\\{p,t)\ 2 = 1, (29) 

and, therefore, the total probability density of finding a particle in the given 
space-time point ( x , t) 


\^(x,t)\ 2 = J dp\i/j(r,t)\ 2 

= \4>>(x,t)\ 2 + |^<(x,t)| 2 + 2Reil> > (x,t)ip* < (x,t), (30) 

as it follows from (l26l) . and the functions ip > ( < ' ) (x,t) are determined by Eqs. 
(El where '0>(<)(£'; p,t) is replaced with ip > ( < ' ) (E) (see (1751) 1. 

A physically relevant situation occurs when the initial wave function vanishes 
at x > 0 (well localized within the x < 0 half-line) because the propagator (l20l) 
transmits this function from the x’ < 0 region to the x > 0 or x < 0 regions. 
This can be achieved if we define the initial wave function as El at x' < 0 and 
set it zero at x' > 0. It can be shown that when the condition 

IDI » 1 < 31 > 

holds (i.e. when the tail of the initial wave packet El is very small near the 
arrival point x = 0), the Fourier transform of the initial wave packet matches 
the Fourier transform of a cutoff Gaussian wave packet, defined as El at x' < 0 
and zero at x' > 0 (see [22]). 

Generally, both the i/j > (r, t ) and t/>< (r, t) components contribute to the prob¬ 
ability density |^>(r, f)| (see (]26]) 1. We can also assume that 


kia 1 , 


(32) 
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which implies that the perpendicular-to-interface momentum dispersion h/2a 
is much smaller than the corresponding characteristic momentum pi = hki , or, 
equivalently, 


H 2 

2 ma 2 


<C E±, 


E± 


2u2 


h 2 k 

2m 


= Ei- 


h 2 ki 2 

2m 


(33) 


i.e., the energy dispersion H 2 /8ma 2 is much smaller than the perpendicular 
component E± of the incident particle energy £) = (h 2 /2m){k 2 + ). Then 

one can see from (12dl) and (051) that in the case when condition m holds, 
the contribution of the backward-moving term ip < (r,t) to the probability den¬ 
sity is significantly smaller than that of the forward-moving term ip > (r,t), and, 
therefore, in the first approximation the former can be neglected. Thus, the 
backward-moving term ip < (r, t ) is not essential in the quasi-classical approxi¬ 
mation when both inequalities (1311) and (1551) are satisfied and, therefore, the 
particle scattering at the potential 0 is associated with the wave packet m 
characterized by a well-defined location relative to the potential and well-defined 
momentum. However, if the inequality (1521) (or (1551) 1 is violated, then both the 
forward- and backward-moving components of the wave function (1241) equally 
contribute to the probability density |^>(r,t)| . In this case the quasi-classical ap¬ 
proximation is not relevant and the particle is associated with the well-localized 
wave packet which has the broad perpendicular-to-interface momentum (energy) 
distribution. 


4 Stationary case and numerical modeling 

We will consider the probability density \i/j(x,t)\ 2 (151)1) . It is convenient to shift to 
dimensionless variables. As seen from (l24l) . there is a natural spatial scale d, an 
energy scale Ed = H 2 /2md 2 (the energy uncertainty due to particle localization 
within a barrier of width d), and a corresponding time scale td = h/Ed . Then, 
using (l24l) and (l28l) (with C\\(p, t) = 1), we can obtain the wave function ip(x, t ) 
in the different spatial regions resulting from the evolution of the initial Gaussian 
wave packet (l27l) in the presence of the potential barrier 0. Thus, the one- 
dimensional wave function, following from (1241) and needed for the calculation 
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of the probability density (1301) , in the dimensionless variables is 
ip(x,t) =ip > (x,t) + ip < (x,t),t> t 0 , 


4>>{x,t) = 


1 a 1 / 2 


,iki. 


oo 


dE 


o -iE(t-t 0 


v / 2(2tt) 3 / 4 d 1 / 2 J E 1 / 4 (E - A) 1 / 4 

Xe -(V&-Vfi ±_) 2 v 2 e W E-A(x-1) e -i\Z~Exi > ^ 


h(E) 


'fp<{x,t) = 


Ip>(x,t) = 


1 a 1 / 2 


AkiXi 


dE 


V / 2(27t) 3 / 4 d 1 / 2 j £!/ 4 [(£- A) 1 / 4 ] 




1 a 1 / 2 




d# 


= -iB(t-t 0 ) E-t/a; 


V / 2(2tt) 3 / 4 d 1 / 2 7 S 1 / 4 (E - C/) 1 / 4 

+r'(£)e- < ^ s ]e-(^-^) 252 e- i ^ I 0 < 2 < 1, 


4><(x,t) = 


1 cr 1 / 2 

v / 2(2tt) 3 / 4 d 1 / 2 ' 




d£ 


-iE(t-t 0 ) \^i^^ p i\/~E-Ux 


J E^KE- U) 1 / 4 } 

+r \E) e - i VE^U^* e -(y/E+\/E ± _) 2 a 2 e iy/Ex ii o < J < ! 

I 1 CT 1 / 2 

1 p>{X,t) = 


„iki: 


f dE ^ -iE(t-tg) ^i'/Ex 


\/2(27r) 3 / 4 d 1 / 2 7 

+r(£)e- iv/ ^ i ]e- ( ^-^ )252 e- , ^,i < 0, 


ip<(x,t) = 


;i/2 


AkiXi 


v / 2(27r) 3 / 4 dV2 J ^ 


oo ^ 

[ dE c -iE(t-to) | e -iV~fix 

J \/f, 


+r*{E)e i ^]e~^ + ^^ 2cy2 e^ Xi ,x < 0, 


(34) 
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where 


m 

t'{E) 

r'(E) 

r{E) 

d{E) 


4 £’ 1 / 4 '\/E — U(E - A) 1 / 4 e iy /~E~u 

W) ’ 

2 E x ' 4 (E - Uy/ 4 {VE- A + VE-U) 

W) ’ 

2 E X I\E - Uy/ 4 {VE-U - Ve - A)e 2i ^~^ 

W) ’ 

(Ve- Ve- u)(Ve -a + Ve-u) 
-(Ve + VE-U)(V'E- A- \/e - U)e 2l 'fz-u 

(>/l + \Je~U)(\Je~A + \Je-U) 

(VI- ^E-U){^E-A~^E-U)e 2i ^~^, 


(35) 


and E =~E/E d , U= U/E^A = AJE d , E ± = E ± /E d , E ± = H 2 k 2 /2m, 
t = t/t d , to = to/t d , a = a/d , x = x/d, Xi = Xi/d. The conditions (l3lT) and (1331) 
read in the dimensionless variables, correspondingly, 


\xi\ » 2a, Ex » 1/a 2 . 


(36) 


It is instructive to consider first the limiting case defined by the second 

inequality (l36l) . In this case, the forward-moving terms ip>(x,t ) in Eqs. (l34l) 

give the main contribution to the total wave function, i.e., i l>(x,t) ~ ip>(x,t). 

Also, the integrals over energy in ^>>(x,i) (|34l) can be asymptotically evaluated 

at A = Eid 2 1 due to the fact that the contribution to these integrals mainly 

comes from the energy region E « E±. In this case, the wave functions ip > (x, t) 

reduce (in the first approximation with - 1 <C 1) to the stationary (for E = 

E j_ <j z 

Ej_) results, oscillating with time as exp[— iE±(t — £ 0 )]. Thus, if we present Eqs. 
(1571) for ip > (x, t) as 


OO 


^>0M) = j (p(x,Xi\E)exp[-iE(t-to)]exp[\f(E)]dE, 
o 

A = £ ± a 2 »l,/(E) = -(VI- \fK) 2 /E ± , (37) 

where ip(x, Xf, E) stands for any integrand in (1571) multiplied by exponentials of 
(1371) . the asymptotic value of (1371) is 

ip>{x,t) ^ ^-J^-\/So_(p(x,Xi-, E±) exp[-iE±(t -to)]. (38) 

a 
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Accordingly, this stationary result leads to the square modulus of the wave 
function |'0 > (:r,f)| , defined by Eqs. (1M1) . which is independent of time. For 
the case of the potential well (U < 0) as well as for the potential barrier (U > 0), 
we obtain at E± > A (in the original non-scaled variables) 


, , , m2 1 16El \E± - U\ ^ J 
W>(x)\ =~^= - , „„ ,, 2 - ,x>d, 


\ 4 ’> 0*01 = 


V^cr \d(E ± )\ 
1 16E ± 


E ± - A - (U - A) cos 2 [y/2m(E± - U)/H 2 (x - d)] 


V^a \d(Ej_)\ 2 

\d{E±)\ 2 = 4 (y/F[ + \/E± - A) 2 (E± -U) + 4U(U - A) sin 2 [ V / 2m(F J _ - U)/h 2 d\ 
E _l > H 2 /2mcr 2 . (39) 


, 0 < x < d, 


Note that when a particle tunnels through a barrier (E± > 0, U >0, E± < U), 
cos[y / 27n(£ , _L — U)/h 2 (x — d)] and sin[y / 2m(Fj_ — U)/H 2 d] in (15^1) should be 
replaced with cosh[^/2m(L/ — Ej_)/h 2 (x — d)] and i smh[\/2m{U — Ej_)/h 2 d\, 
respectively. 

Formulae (l39l) provide the spatial dependence of the wave function square 
modulus at different spatial regions relative to the potential area for the station¬ 
ary case, when the initial wave packet m is characterized by an extra narrow 
distribution in the energy (perpendicular-to-interface momentum) space. Thus, 
in this approximation, the transmitted probability density (x > d ) is constant 
in space, while in the potential region (0 < x < d) we have the oscillating 
interference pattern (for E± > U). 

The picture before the potential (x < 0) is more complicated and results from 
the interference of the incoming and reflected waves. The corresponding formula 
becomes simplified for the resonant case, when yj2m(E± — U)/h 2 d = ^m (n is 
the integer), and is given by (E± > U, E± > A) 

l ’ M “° |2 = Wm. + m - - A S in 2 (V2mE 7Jh^)U < 0. 

Ej_ H 2 /2ma 2 , \J2m(E±_ — U)/h 2 d = 7m. (40) 

The oscillating interference picture given by (1401) is caused by the earlier-mentioned 
fact that in the case of an asymmetric potential (A ^ 0), the reflection amplitude 
r(E) ^ 0 for the resonant energies E (see (fT7l) ~). From Eqs. (l39l) and (1401) we 
see that the norm |'i/>>(0)| = (^e ± +^e±-A)' 2 P°t en ti a l left bound¬ 

ary x = 0 is transmitted at the resonance condition y/2m(E± — U)/h 2 d = nn 
to the region x > d beyond the potential. Only for a symmetric rectangular 
potential (A = 0) the reflection amplitude r(E) = 0 for the resonant energies 
and there is only the probability density |^>(a;)| = l/\/27r<7 (x < 0) stemming 
from an incoming wave and arriving to the x > d area. Thus, the dependence of 
the constant in space transmitted probability density (1391) versus the potential 
amplitude U will exhibit the oscillating (at E± > U ) pattern beyond the barrier 
(x > d) with an amplitude which is greater for the asymmetric barrier (A ^ 0) 
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as compared to the symmetric one (A = 0). The same is true for the oscillating 
^-dependence of |^>>(x)| inside the potential area (0 < x < d). 

The time dependence of the probability density \ip(x, t) f exhibits itself only 
when there is a sufficient momentum dispersion, as follows from Eqs. (1M1) . On 
the other hand, a sufficient momentum dispersion, when E±d 2 ~ 1, leads to a 
nonnegligible counterintuitive contribution of the backward-moving components 
of the wave packet to |^(a;,t)| 2 . The spacetime evolution of the scattering 
process can be visualized by numerical evaluation of the probability density 
(l34l) of finding the particle in the scaled space-time point (x,t). We 
will focus on the influence of the wave packet backward-moving components and 
the potential asymmetry parameter A on the particle dynamics. As mentioned 
earlier, the asymmetric rectangular potential can model the potential profile of 
the magnetic threelayer when it is switched from the parallel configuration of the 
magnetic layer (modelled by the symmetric potential profile with A = 0) to the 
antiparallel orientation. For the case under consideration, when the particle, 
associated with the Gaussian wave packet, moves towards the potential m 
from the left, one can expect that the influence of the asymmetry parameter A 
(defining the height of the right potential step of (fljl) will be more pronounced 
if the contribution of the backward-moving components of the wave packet is 
essential (the numerical evaluation confirms this expectation). 

To make the dynamics of the wave packet more particle-like, we accept the 
condition of the narrow wave packet, a < 1, and put to = 0. For an electron and 
the potential width d = 10 -7 c?n (lnm), the characteristic energy Ed ~ 3-10~ 2 ev 
and the characteristic time id ~ 2 • 10” 14 s. In accordance with the accepted 
conditions, we will posit E± = 10 2 , Xi = —10, and d = 1/3 or d = 0.1. We 
choose C7 = 10 in the case of a potential barrier (over-barrier transmission), and 
U = —10 2 for a potential well. We will compare two cases: d = 1/3, when the 
second inequality (l36l) is satisfied and the backward-moving positive energies 
components of the initial wave packet are not essential, and d = 0.1, when their 
contribution matters. The dimensionless time interval t = 0.1 -r 1.5 is chosen 
from a simple estimation for the average scaled time tj /td that it takes a particle 
with the initial energy E± = 10 2 to reach the potential starting from the point 

Xi = -10: ti/t d = \xi \ m/hkitd = \xi\ /2 \f^± = 1/2. 

Figure 1 shows the probability density k/(a;,f) of finding the particle at 
x = 1, i.e. on the right-hand side of the barrier ([]]) (U > 0), as a function of t 
and A changing from A = 0 to A = E±/ 2 when a = 1/3. Figure 2 shows the 
same function for a = 0.1. We see that in the case when the contribution of 

the backward-moving components of the wave packet is important ( d = 0.1), 

~ . 2 

the time distribution of finding the particle beyond the barrier ^>(l,t) for the 
asymmetric potential is essentially different from that for the symmetric one: 
Beginning from the value of the asymmetry parameter A ss 20, this distribution 
becomes more broad and pronouncedly nonmonotonic for A > 20. 

For the case of a potential well with U = —10 2 , we numerically evaluated 
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_ 2 

Figure 1: Probability density distribution |i/(l,f)| on the right-hand side of 
the barrier as a function of time and asymmetry parameter A for the narrow 
energy distribution of the initial wave packet (a = 1/3). 



f 


Figure 2: Probability density [^(l, ^) | 2 as a function of t and A for the broad 
energy distribution of the initial wave packet (cf = 0.1). 
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Figure 3: Probability density |^>(a;,f)| inside the symmetric well (A = 0) for 
the broad energy distribution of the initial wave packet (cr = 0.1). 


inside the well (x = 0 4-1) for the case of the broad energy distribu¬ 
tion of the initial wave packet (a = 0.1) and the asymmetry parameter A = 0 


and A = E±/2 


U 


72- 


Figure 3 shows the interference pattern inside the 


symmetric well which differs sufficiently from the stationary square cosine type 
picture, given by Eq. (l39l) (for A = 0). It is seen that the amplitude of this 
pattern grows with time from zero to the maximum value (reached approxi¬ 
mately at t = 0.5) and then again diminishes to zero, thereby showing the finite 
time during which a particle exists in the well region before leaving it either 
for the region before (x < 0) or beyond (x > 1) the well. We also see that the 
interference pattern of t)\ is more structured in space and time. These 
changes in the probability density distribution result from the influence of the 
backward-moving components of the wave function ip < {x,t), which is essential 
for the considered case of sufficient energy dispersion (E±d 2 = 1). In Fig. 4, 

we see the influence of the asymmetry parameter (A = 50) on that probability 
. __ . 2 

density ^(x, t) inside the asymmetric well. The calculated distribution ex¬ 
hibits a very structured and pronouncedly nonmonotonic interference pattern 
in space and time compared with that displayed in Fig. 3. 


5 Dwell time 

For a finite spatial interval, the so-called dwell time, i.e. the average time spent 
in this interval by a particle described by the packet ip(x,t), is customarily 
used. The dwell time in the potential region (0, d) can be defined in the three- 
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1.0 


Figure 4: Probability density |i/>(a;,f)| inside the asymmetric well (A = 50) for 
a = 0.1. 


dimensional case for the potential m as 


T (0, d) 


lim 

t —to—too 


dt dx / dp \^(x, pd)Y 


to 


(41) 


Substituting the third and fourth lines of (El) into the definition (TUI) , we obtain 
for the initial wave packet m 


(E) 


(42) 

where the functions ip > r < \(E) are defined by C51) and the equality (El) is taken 
into account. 

We see that, again, the dwell time is determined by the forward- and backward- 
moving components of the initial wave packet as well as their interference. The 
entire range of energy (0 4 - oo) contributes to the dwell time. It is not diffi¬ 
cult to get from (1421) . (1171) and m that the forward- and backward-moving 
components of the dwell time are 

OO 

t>«)M= JdET H<) (E-d). (43) 

0 


r(0,d) = r>(0 ,d) +r<(0,d) + 2Rer> 

oo d 

(0, d)= dE / dx 


J J 

0 0 

V v u{ E) 


oo d 


2 

r >< (0,d) = / dE / dx 

v{x,E) 

i>>(i 

J J 

0 0 

_y/v u (E)_ 

<p(x, E) = t'{E)e lk " x + r\E)e~ ik - x , 


OO 

d 

[ dE 

/ dx 

J 

0 

J 

0 


ip(x, E) 


VME) 


< 
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The per unit energy interval energy-dependent dwell time r>(<) ( E ; d) caused by 
the forward- (backward-) moving component of the initial wave packet ^>(<) (E) 
is for E > A (/ca = i/2 m(E — A)/h 2 is real) 


r>«)(E-,d) = t(E; d)\ip >{<] (E)\ 


t(E;d) = 


k 2 k u d(k 2 + — Kq sin(2 k u d) 

v u (k + k^) 2 k 2 + k^Kg sin 2 (k u d) 


,E > A > 0, 


(44) 


where K 2 = k\- k 2 u = ^-(U - A), k% = k 2 - k 2 u = U . For E < A, A > 0 


(/ca = ik&, k'A = y / 2?n(A - E)/h 2 ) 


k 2 k u d(k 2 + &a) — Kq sin(2 k u d) + 4/c u A’a sin 2 (k u d) 
Vu (k + k,A.) 2 k 2 + IcqA'q sm 2 (k u d) + k'^kuk^ sin(2 k u d) 


t(E;d) = 


,E < A, A > 0, 

_ _ (45) 
where I\ \ = /ca — k 2 = jjr{U + A — 2 E). Both Eq. (El) and Eq. El) are valid 

for k u real as well as for the imaginary k u = ik u , k u = y/2m(U — E)/h 2 when 
E < U (barrier). Note that ^ > ( < )(£’) has the dimensionality of the inverse 
square root of energy (see (1281) 1. and thus the expression for t(E;d) has the di¬ 
mensionality of time and represents the generalization of the energy-dependent 
dwell time obtained earlier by Buttiker [23] to the case of the asymmetric rect¬ 
angular potential ([TJ (if A = 0, Eq. El I reduces to the Buttiker result). 

The per unit energy interval interference dwell time which follows from 
can be written as 


t ><: (E ; d) = Re ■ 


k ( k 2 + k\) sin(2 k u d) + 2(/c 2 — k\)k u d — Aik u kA sin 2 {k u d) 


[i(kkA + k 2 ) sin (k u d) — k u (k + A’a) cos (k u d)] z 


V>>(£)V><(£) 


Rer><(0,d) = /dEr><(E;d). 


(46) 


Note that Eq. El) holds for both k u real {E > XJ) and imaginary k u = ik u 
(E < U), as well as for both the real A;a (E > A) and imaginary A’a = //ca 
(E < A). 

We see that the total per unit energy dwell time 


T(E-d)=t(E-d) |^>(e ;)| 2 + |<me )| 2 


■ 2r >< (E;d),r(0,d) = J r{E\d)dE 
0 


(47) 

is generally defined by both the forward- and backward-moving components of 
the initial wave packet as well as their interference. For the resonance energies 
satisfying the condition k u d = nn (n is integer, k u is real), taking place in the 
cases of U < 0 and U > 0 (when E > U), Eqs. El) - d47l) reduce, e.g. for /ca 
real (E > A), to 


t(E; d) = 


d ( 2 k 2 (k 2 + k%) 


v ( kl(k + /ca ) 2 


\^>{E)\ 2 + \^<{E)Y 


4 Re 


k 2 (k 2 u -k 2 A ) 
k 2 (k + /ca ) 2 


[4’>(E)ip* < (-E 1 )] 
(48) 
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where d/v{E ) is the time that it takes for a particle with the energy E to 
propagate through the spatial range d in the absence of a potential. Thus, the 
expression in the curly brackets in (1481) shows the difference between the dwell 
time in the range of the potential and the ’’free” dwell time d/v(E). 

From the above it follows that, generally, the dwell time depends on the 
energy spectrum of the initial wave packet ip > ( < \(E) and cannot be realistically 
defined, e.g., simply by t(E\d) (1441) or (l45l) . Further, we will use 4’ > { < ){E) (|28|) . 
defined for the Gaussian initial wave packet, and shift to the dimensionless 
variables defined in the previous section. As a result, we obtain from Eqs. dm > 

- (S3 


oo 

r(0 ,d) = jdEr(E\ d),r(E; d) = T>(E;d) + r < (E;d) +?><(£(; d), 
o 

~ ,X. ^ ~ ^ td* 1 2 VE - U(2E-U-A)-(U-A) sin(2 \/E - U) 

t >\E, d) + t<[E, d) _— — — -— ~ ~ ~ ~ ~ rsss -— 

2V2t r yjE-U (v E + V E - A ) 2 (E - U) + U{U - A) sin 2 (V£; - U) 

x |exp[-2(/iI - \fEfo 1 } + exp[-2(-/§T + \/l) 2 d 2 ]|, E > A > 0, 

T><{E ' ,d) = v^ x 

1 (2 E-U- A) sin(2 \JE -U) + 2\Je~ U{ A -U)-li\/E- A y/E - U sin : 

VE-U [i(VEy/E -A + E-U) sin (\/E - U) - (Ve + V E - A)\/e — U cos(V 

x exp[— (\JE± — \/e) 2 i t 2 ] exp[— (\JE± + \/e) 2 (j 2 ], (49) 


where the characteristic time td /2 = h/2Ed = \/~E = v ^ ± -j \J E±, i.e. it is 

the time spent in the region of the potential width d by a ’’free” particle with 
the energy E = Ed (E = 1), and thus t(E\ d) has the dimensionality of time (for 
brevity, we do not show Eq. (l45l) in the dimensionless variables). The relative 
contribution of the forward- (backward-) moving components r>(<)(Fl;d) and 
interference term r><(E; d) to the dwell time r(0, d) (l49l) depends on the value 
of the parameter Eid 2 . If the second inequality (l36l) is satisfied, i.e., E±d 2 1, 
the contribution of the backward-moving and interference terms to the dwell 
time (1491) is much smaller than that of the forward-moving term r>(E; d), and, 
therefore, the former terms may be ignored in the first approximation in the limit 
given by m- Moreover, the integral of r> {E-. d) over E can be asymptotically 
estimated due to the sharp maximum of the integrand at E = E±. The result 
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IS 


r(0,d) ~ r>(0,d) = 


E± > A > 0, 


E± 


2 a/ E± - U(2E± — U — A) — (U — A) sin(2 JE ± - U ) 


>{E±) 


E± — U (\JE± + yjE ± - A) 2 (E± -U) + U(U - A) sin 2 (^-Ex - U) 

(50) 


which coincides with Eq. (l44l) for t(E±\d) written in the dimensionless vari¬ 
ables. It should be stressed that this result represents only the first term of 

OO 

the asymptotic expansion of J dEr > (E; d) with a small value of the parame- 

o 

ter 1/E±d 2 , i.e. for an initial wave packet characterized by an extra narrow 
momentum distribution. 

For the resonance energies satisfying the condition k u d = nn ( n is integer, 
n ^ 0, k u is real), which reads in the dimensionless variables as E± — U = n 2 n 2 , 
the relative to the ’’free” dwell time d/v(E±) expression (T5U1) reduces to 


T r > {Q,d)v(E 1 _)/d = 


2E ± (2Ej - U - A) 


2E±(Ej — A + 7r 2 n 2 ) 


E ± + JE ± -A) 2 (E ± -U) 


At U < 0 (dwell), the inequality 


U 


Ej_ + y E _l — A) 2 TT 2 n 2 

(5i) 

= 7r 2 n 2 — E_i > 0 should be satisfied 


[n is bottom-limited), and when 7r 2 n 2 3>> E± > E± — A > 0, which is the 


case for large enough 


, the asymptotic relative resonant dwell time m 


u , 

approaches 2Ej_/(\JE± + \J E± — A) 2 . This value is greater than 1/2, to which 
the values of the high order resonances of the dwell time reduce for a symmetric 
potential (A = 0). Thus, the greater the asymmetry parameter A, the greater 
the amplitudes of the dwell time resonances. For U > 0 (barrier), the condition 
U = E± — ir 2 n 2 > 0 should hold (E± > U, n is restricted to the small values 
defined by E±), and at Tr 2 n 2 <C E± — A < E± the dwell time (1511) behaves as 

2E±(E± - A)/{yfl[^+yjEj_ - A) 2 t r 2 n 2 . 

If E± — U —> 0 (reverse points in classical physics), which can happen only 
at U > 0, the asymptotic relative dwell time r>(0, d)v(E±)/d (l50l) reduces to 

4 E± AE± 


r > (0,d)v(E ± )/d = 


> 


E ± + JE ± -A) 2 + U(U-A) 


4E± + U 2 


< 1, E±—U —> 0, 


where the value AE±/(AE± + U 2 ) corresponds to a symmetric potential, 


(52) 

i.e., 


the dwell time (1521) in the asymmetric case is larger. In particular, at U = 0 
this dwell time for A > 0 is larger than the ’’free” dwell time d/v{E±) in the 
absence of a potential (U = 0, A = 0). 
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r>(0,d)v(£J/d 



Figure 5: Dependence of the asymptotic dwell time (1501) on the well/barrier 
symmetric potential (A = 0). 


It is interesting to plot the dependence of the relative dwell time (15(1 on 
the amplitude of the potential U, which changes from negative values (the well) 
to positive ones (the barrier), for a symmetric (A = 0) and an asymmetric 
potential. Formula (1501) is valid for both the U < 0 and U > 0 cases, and in the 
latter case, when E± — U < 0, Eq. (EHD transfers to 


T > (0,d)v(E±)/d 


E± 


2\JU — E ± (U + A - 2E ± ) + {U - A) sinh(2 sJU- E L ) 


U - E ± (y/E .l + yjE _ L - A ) 2 (U - E _ L ) + U(U - A) sinl - E _ L ) 


(53) 

Fig. [5] shows the U - dependence of r>(0, d)v(E±)/d (l50l) for a symmetric 
potential (A = 0) in the broad range of U = —2 x 10 3 -f- 2 x 10 2 and Ej_ = 
10 2 . One can see the series of resonances at U < 0, the amplitudes of which 
approach 1/2 for big enough n, n 2 E±/ir 2 , while at U > 0 there is a limited 
series of resonances with n 2 < E±/tt 2 (for E± > U) with the larger amplitudes 
because a particle moves more slowly in the presence of a potential barrier 
than in the region of a potential well. Fig. [G] shows the same U - dependence of 
r>(0, d)v(E±)/d (l5Ul) for A = 90. We see an essential increase of the resonances’ 
amplitudes inside the well and other details which display the influence of the 
potential asymmetry on the dwell time in correspondence with the analysis given 
above. 


26 












r > (0,d)v(EJ/d 



Figure 6: Dependence of the asymptotic dwell time (1501) on the well/barrier 
asymmetric potential (A = 90). 


6 Summary 

We have applied the MST to the calculation of the propagator which exactly 
resolves the time-dependent Schrodinger equation for a particle in the pres¬ 
ence of a one-dimensional rectangular asymmetric well/barrier potential (fT|). 
This approach, based on the obtained effective potentials 0, ©, which are 
responsible for reflection from and transmission through the potential steps, is 
alternative to the matching procedure conventionally used for solving the sta¬ 
tionary Schrodinger equation. The advantages of this MST approach are: A 
natural picture of the considered processes in terms of a particle scattering at 
the potential jumps (in contrast to the traditional wave point of view); The 
time-dependent picture of the quantum effects of particle reflection from a po¬ 
tential well and particle transmission through a potential barrier; The natural 
decomposition of the Schrodinger equation solution into the sum of the forward- 
and backward-moving terms (with no use of the evanescent states [3]), which 
takes into account that the initial wave packet, confined to a restricted spa¬ 
tial area and representing a particle moving towards a potential, contains both 
the positive and negative momentum components. Aside from being related to 
the fundamental issues of quantum mechanics, the obtained results can be also 
important for the kinetic theory of nanostructures, where the considered rect¬ 
angular potential 0 is often used to model the potential profile in the magnetic 
nanostructures utilized, e.g., in spintronics devices. 

The obtained probability density | t) | 2 of finding a particle in the space- 
time point ( x , t), when it initially was located in some spatial region and moved 
in some direction, is generally defined by the probability density correspond¬ 
ing to the wave component moving in this direction \if> > (x,t)\ 2 as well as by 
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the probability densities related to the backward-moving component \ip < (x, t)\ 2 
and the interference of both 2R e['i/j > (x,t)il)<(x,t)]. For the case of the initial 
Gaussian wave packet, we have shown that the contribution of the backward- 
moving component to the probability density \tp(x, t)| is small when the initial 
packet is characterized by a narrow energy (momentum) distribution, which is 
characteristic of the quasi-classical approximation for a transport phenomenon. 
We calculated, in this case, the asymptotic time-independent values of |^>(a;)| 2 
in the different spatial regions relative to the potential area. This situation (ex¬ 
tra narrow energy distribution) actually corresponds to the stationary case with 
no energy dispersion. Thus, the transmission through and reflection from the 
potential well/barrier can be described as a function of time only when the mo¬ 
mentum (energy) dispersion of the initial wave packet is significant (accordingly, 
the wave packet spatial localization is narrow). But in this case, the counterin¬ 
tuitive (non-classical) contribution of the backward-moving components of the 
wave packet should be accounted for. This rather paradoxical quantum me¬ 
chanical result reveals itself in the problems connected to measuring time in 
quantum mechanical effects. 

Using the exact result for \i/)(x, i)| 2 , we have numerically plotted the time 
distribution of finding the particle beyond the barrier (U > 0), |^(l,t)| and 
found that, when the contribution of the backward-moving wave packet com¬ 
ponents is important (broad wave packet energy distribution), the influence of 
the potential asymmetry can be essential (Figs. 1,2). Plotting \ip(x, t)\ 2 in the 
well (U < 0) region, we showed that the backward-moving components of the 
wave packet fundamentally change the probability density, when the initial wave 
packet is broad enough in the energy (momentum) space, and the asymmetry 
of the potential well adds more to the structure of this spacetime distribution 
(Figs. 3,4). 

The obtained solution is applied to the calculation of the particle time dwell 
time within the potential area. Again, the forward- and backward-moving com¬ 
ponents of the obtained exact wave function contribute to the particle dwell 
time. For a narrow momentum distribution of the initial wave packet, the an¬ 
alytical asymptotic value of the main (in this case) term contributing to the 
dwell time in the potential region, caused by the forward-moving probability 
density \ip > {x,t)\ 2 , was obtained and plotted as a function of the potential am¬ 
plitude U changing from the negative (well) to the positive (barrier) values. 
The series resonances displayed in Figs. 5,6 show the essential influence of the 
potential asymmetry on the particle dwell time. These results generalize the 
known Buttiker results [23] for the dwell time. 
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